Draft version May 11, 2010 

Preprint typeset using WT^^ style eniulatcapj v. 11/10/09 



COMPACT SUNYAEV-ZELDOVICH 'HOLE' IN THE BULLET CLUSTER 

S. S. Malu^ and R. Subrahmanyan 

Raman Research Institute, Sadashivanagar, Bangalore 560 080, India 



o 

(N 

>>■ 



o 
u 

6 



en 



in 

o 
o 



X 



M. H. Wieringa 

Australia Telescope National Facility, CSIRO Astronomy and Space Science, Locked Bag 194, Narrabri NSW 2390, Australia 

AND 

D. Narasimha 

Tata Institute of Fundamental Research, Mumbai 400 005, India 

Draft version May 11, 2010 

ABSTRACT 

We present 18 GHz observations of the Bullet cluster from the Australia Telescope Compact Array; 
in particular, a high angular resolution measurement of the substructure in Sunyaev-Zel'dovich Effect 
(SZE). We report the discovery of a compact SZE 'hole' in the galaxy cluster. The SZE hole does not 
correspond to any bright feature in X-ray, optical or lensing maps; in general, the relatively deeper 
SZE features appear to avoid the regions with the most intense X-ray emission. These imply that the 
gas pressure distribution differs significantly from the distributions in gas emission measure, galaxy 
and dark matter distributions. This has implications for the gas physics and evolution in the cluster 
merger event. SZE displaced from X-ray centres implies that modeling cluster dynamics is non-trivial; 
our observations indicate that our current lack of understanding cluster merger astrophysics may be 
a limitation to modeling the cosmological distribution in SZE cluster counts and the cluster SZE 
contribution to small-angle cosmic microwave background (CMB) anisotropy. The SZE distribution 
in the western parts of the cluster are cospatial with the radio halo indicative of a common origin for 
the hot and relativistic electrons in the turbulent wake of the Bullet. 
Subject headings: cosmic background radiation — galaxies: clusters: individual (IE 0657-56, RX 

J0658-5557) — intergalactic medium — radio continuum: general — techniques: 

interferometric 



1. INTRODUCTION 

1E0657— 56, known as the 'Bullet Cluster', is one 
of the hottest known clusters that has been well- 
studied over the last decade for a variety of reasons; 
name ly, the existence of a strong radio halo (iLiang et al' 
20001 ) ■ the Sunyaev-Zel'dovich effect ([Halverson et al. 



and references therein), its X-ray brightness 



( Markevitch et al.]|2002|) , though most notably in provid 
ing the most direct proof of the existence of dark matter 
(jClowe et al.ll2006l ). It is a cluster collision/merger event 
at z ~0.296, with the larger, westward cluster being ^10 
times the mass of the smaller 'bullet'. 

Observation s fro m the South Pole Telescope 
(jLueker et al.l I2009D show that on arcminute and 
smaller scales observed CMB anisotropy power is signifi- 
cantly less than levels expected from models of SZE from 
a cosmological distribution of clusters. This discrepancy 
may be an indication of our lack of understanding of the 
environmental impact of cluster mergers. This paper 
presents sub-arcmin angular resolution observations of 
the SZE in the Bullet cluster, which demonstrates that 
the relation between SZE and X-ray emission is complex 
and requires additional physics that may be related to 
the impact of the merger on intracluster gas. 

2. RADIO OBSERVATIONS AND IMAGING 
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The Australia Telescope Compact Array (ATCA) is 
a radio interferometer with six 22-m antennas, five of 
which may be positioned on stations along a 'T'-shaped 
rail track that is 3-km along E-W and 214-m along N-S. 
Observations of the Bullet cluster were made in a pair 
of 2-GHz bands: a '17-GHz band' covering frequencies 
16-18 GHz and a '19- GHz band' covering the range 18- 
20 GHz. Each of the 2-GHz wide bands were subdivided 
into 2048 frequency channels. All observations were in 
full polarization mode and recorded multi-channel con- 
tinuum visibilities. 

X-ray images of the cluster show two peaks separated 
by ^ 1.'5 and the detected extent of X-ray emission is 
~ 6' (jAndreani et al.HlQOQt ) suggesting that hot gas per- 
vades a significant area between and beyond the con- 
stituent clusters. Since the FWHM field-of-view of the 
ATCA antennas is 2f6 at 18 GHz, we mosaic imaged the 
Bullet cluster using pointing centres at (J2000 epoch co- 
ordinates) RA: 06'^58'"20^ DEC: -55°56'28" and RA: 
06h58'"30^ DEC: -55°56'28". These two pointing cen- 
tres correspond to the X-ray brightness peaks and are 
separated by 0.54 of the FWHM of the ATCA antenna 
primary beams at 18 GHz. 

Observations were made in each of two ATCA array 
configurations: H168 and H75 that have baselines up to 
168 and 75 m respectively; the former to enable subtrac- 
tion of unresolved continuum sources, the latter to pro- 
vide enhanced surface brightness sensitivity for imaging 
of the SZE. A journal of the observations is in Table [TJ 



Table 1 






Journal of the ATCA observations. 


Array configuration Observing time 


(flours) 


Date 


H168 12 




2009 April 25 


H168 12 




2009 April 26 


H75 8 




2009 June 28 


H75 8 




2009 June 29 


H75 8 




2009 June 30 



In each of the observing sessions, antenna pointing cor- 
rections were updated every hour using a 5-point offset 
pattern observation on PKS B0537— 441, the unresolved 
cahbrator PKS B0742— 56 was observed every 10 min to 
monitor and correct for amphtude and phase drifts in the 
interferometer arms, and PKS B1934— 638 was observed 
once every session as a primary cahbrator to set the ab- 
solute flux density scale. During the 10 mins between 
successive calibrations, the two pointing positions were 
alternately observed for 80 sec each. Visibilities were 
recorded with 10 sec averaging. 



The interferometer visibilities were examined, calibrated 
and imaged using MIRIAD; all image processing were 
also accomplished using utilities in this software pack- 
age. The visibility data in the 17 and 19 GHz bands 
were separately edited for interference and calibrated be- 
fore bandwidth synthesis imaging. Visibility data in each 
of the 2 GHz wide bands were recorded over 2048 fre- 
quency channels, and 50 channels at each of the band 
edges were excluded from analysis to avoid data in fre- 
quency domains where signal path gains are relatively 
low. Frequency channels that appeared to have relatively 
large fluctuations in visibility amplitude owing to hard- 
ware faults in the digital correlator were also rejected 
prior to calibration and imaging. 

2.1. Calibration 

Adopted fluxes for the primary calibrator 
PKS B1934-638 were 1.146 and 0.992 Jy in the 
17 and 19 GHz bands respectively; the spectral index 
a (defined as S^, oc v", where Si, is the flux density 
at frequency v) was adopted to be —1.33 in both 
bands. Outliers in the amplitudes of visibility data 
on PKS B1934— 638 were rejected — removing 15% of 
data — and the reliable visibilities were used to set the 
absolute flux density scale as well as determine the in- 
strument bandpass calibration. When calibrated for the 
bandpass, the visibility amplitudes of PKS B0742— 56 
showed continuity across the 4-GHz observing frequency 
range and a trend consistent with a single power-law: 
this was a check of the bandpass calibration. Drifts of up 
to 30° were observed in the interferometer arms over the 
observing sessions: the calibrations for the time-varying 
complex gains in the antenna signal paths as well as 
the calibrations for polarization leakages were derived 
from the visibilities on PKS B0742— 56. The rms phase 
variations in the antenna signal paths within the 1-min 
calibrator scans was within 10°, indicating that short 
timcscale atmospheric and instrumental phase cycling 
would not result in amplitude attenuation of more than 
3%. 



2.2. Imaging 

Since no circular polarization is expected, Stokes V is 
expected to be consistent with thermal noise. Therefore, 
at times and frequency channels where Stokes V visi- 
bilities deviate more than four times rms thermal noise 
in the calibrated visibilities acquired towards the Bul- 
let cluster pointings, data in all Stokes parameters were 
rejected. This Stokes V based clipping was aimed at au- 
tomated rejection of self-generated low level interference. 

The multi-channel continuum visibilities in both 17 
and 19 GHz bands and in both pointings were imaged 
together to form a wide-field multi-frequency synthesis 
mosaic image of the Bullet cluster shown in Fig. [1] (up- 
per panel). The image displayed has been truncated at 
where the attenuation in the mosaic image owing to the 
primary beam — and consequently the response to source 
fiux density — falls below 10%. The mosaic image was 
made with a beam of FWHM 20'.'8 x 14'.'5 at 89° position 
angle (P. A.), and has an rms noise that is 7.4 /iJy beam~^ 
at the centre and increases by factor five to the edges of 
the field of view. The rms noise in the image was es- 
timated from Stokes V data and is consistent with ex- 
pectations based on the system temperature during the 
observations, bandwidth and integration time. The im- 
age displayed has been tapered to make the rms noise 
uniform over the field of view. 

Discrete sources detected by iLiang et"an (|2000( ) at 
1.344 GHz are marked using symbols 'A' through 'H'. 
All these eight sources are detected in the 18 GHz mo- 
saic image at a level exceeding two standard deviation 
in image thermal noise. The brightest source in our 18- 
GHz image is at RA: 06'^58'"37?9, DEC: -55°57'25" and 
is marked 'A' in the figure; this is a LIRG high-z galaxy 
(jRex et al.ll2009( ) and is detected with a fiux density of 
1.69 ± 0.04 mJy at 18 GHz. 

3. IMAGING THE SZE IN THE BULLET CLUSTER 

Radio images of the bullet cluster field inclu de discrete 
radio sources, a radio halo (iLiang et al.ll2000l ). and SZE 
(jHalverson et al.ll2009D that is expected to appear as a 
decrement at the ATCA observing frequencies. For the 
characterization of the SZE decrement, it is essential to 
accurately subtract all discrete sources including the ra- 
dio halo. 

3.1. Subtraction of unresolved continuum sources 

The image in Fig. [1] (upper panel) was made with visi- 
bility weighting close to natural in order to have low rms 
noise and good surface brightness sensitivity. To charac- 
terize first the unresolved continuum sources in the field, 
we constructed a separate image with uniform weight- 
ing of the visibilities that has a factor of two greater 
rms noise but, importantly, de-emphasizes extended halo 
emission and SZE. Sky regions that have image intensity 
exceeding four times the rms noise in Fig. [1] (upper panel) 
were examined for unresolved sources. Apart from the 
sources 'A' through 'H', three sources were detected with 
fiux density exceeding four times the rms noise in the 
uniformly weighted image; these are marked '1' through 
'3' in Fig. [T] (upper panel). Of these, sources '1' and '2' 
are undetected in the 1.344 GHz image and are likely to 
have relatively fiatter spectra. The flux densities of the 
11 unresolved continuum sources in the field were esti- 



mated from the image made with uniform, weighting and 
subtracted from the visibihty data. 

A mosaic image of the Bullet cluster made with close 
to natural weighting and using visibilities from which un- 
resolved continuum sources were subtracted is in Fig. [1] 
(bottom panel) . The beam FWHM and rms noise distri- 
bution and taper are same as for Fig. [1] (upper panel). 




Table 2 

Unresolved continuum radio sources detected in the Bullet 
cluster field 



Source 



RA (J2000) DEC 



Flux Density 



■^^is GHz (niJy) 



A 


06 


58 


38.0 


-55 


57 


25 


1.69 ±0.04 


B 


06 


58 


33.8 


-55 


57 


58 


0.06 ±0.04 


C 


06 


58 


41.4 


-55 


58 


41 


1.4 ±0.4 


D 


06 


58 


23.8 


-55 


56 


41 


0.07 ±0.02 


E 


06 


58 


27.5 


-55 


56 


06 


0.05 ±0.02 


F 


06 


58 


24.5 


-55 


55 


13 


0.11 ±0.03 


G 


06 


58 


19.3 


-55 


58 


41 


0.26 ±0.12 


H 


06 


58 


16.1 


-55 


58 


24 


0.13 ±0.08 


1 


06 


58 


44.6 


-55 


55 


24 


0.5 ±0.2 


2 


06 


58 


36.7 


-55 


54 


57 


0.18 ±0.06 


3 


06 


58 


14.6 


-55 


54 


23 


0.53 ±0.13 



6^58'°50'' 40'' 



lO** 58" 



NOTE S- Letters indicate discrete sources listed inT/iang et al] 
1120001 ) and have the same letters as in their Table 2; numbers 
indicate additional discrete sources from our 18-GHz obser- 
vations. 

06'^58™31!1, DEC: -55°56'27"; this is 18-GHz emission 
from the radio hal o, which was previously discovered by 
iLianget al.) (|2000t) . 
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Figure 1. Upper Panel: 18 GHz mosaic image of the bullet 
cluster. In this representation as well as subsequent images dis- 
played herein the image has been tapered to make the rms noise 
uniform over the field of view; this results in an attenuation of 
unresolved sources by factor 0.02 at the edge. Image rms noise is 
7.4 AtJy beam-l; the beam FWHM is 20'.'8 X 14^5 at 89° P.A. and 
is shown using a filled ellipse in the lower right corner. Contour 
levels are at (—4, —2, 2, 4, 8, 16, 32, 6 4) times the r ms no ise. Lo- 
cations of discrete sources detected by ILiang et al.l 1)20001 ) and in 
this image are marked with cross symbols. 

Bottom Panel: 18 GHz mosaic image of the bullet cluster with 
unresolved continuum sources subtracted. The beam FWHM and 
rms noise are same as for upper panel figure. Contour levels are 
at (—4, —2, 2, 4, 8, 16, 32, 64) times the rms noise. Locations of 
discrete sources detected bv lLiang et al.l 120001 ) and in our 18-GHz 
data are marked with cross symbols. 



The positions of all the sources subtracted have been 
marked in Fig. [1] and in Table [2] we list the positions and 
18-GHz flux densities. A negative feature appears close 
to the location where source 'F' was subtracted; this is 
an extended negative feature that is absent in the higher 
resolution images made with uniform visibility weighting. 
The brightest positive feature in this residual image is 
extended emission with a peak of 40 /xJy beam~^ at RA: 



3.2. Subtraction of the radio halo 

Previous observat ions of the radio halo in the Bul- 
let cluster (jLiang et al. 2000) provide spectral index and 
intensity distribution with sufficient angular resolution, 
and we model the radio halo contribution in our 18 GHz 
image based on their data. There was no evidence 
for spatial variations in spectral index, and the data 
was best fit with a single power-law spectral slope of 
a = -1.3 ± 0.1 between 0.8-8 GHz. 

The 1.344 GHz image of the radio halcH was corrected 
for the telescope primary beam attenuation at their ob- 
serving frequency, an image of the radio halo component 
was isolated, scaled to 18 GHz adopting a spectral index 
of a = —1.3, and attenuated by the mosaic gain func- 
tion at our observing frequency. In Fig. [5] we show the 
18-GHz image of the Bullet cluster using grey scales, he- 
fore any subtraction of the radio halo component, with 
this processed image that represents the expected halo 
component overlaid as contours. 

The 18-GHz image in Fig. [21 which has been smoothed 
to a beam of FWHM 30", clearly shows negative inten- 
sity values towards the central regions of the Bullet clus- 
ter. These SZE features appear to be confined within, 
and bounded by the lowest intensity contours of the halo 
component and are obvious in the 18-GHz image prior 
to any subtraction of the halo component, which would 
only enhance these SZE features. Before subtraction of 
any halo contribution, the 18-GHz image shows a peak 
in flux density at the peak of the halo and positive 18- 
GHz emission towards the eastern parts of the halo. Most 
noteworthy is the negative peak, representing a deep SZE 
'hole', to the south of the peak in the expected halo con- 
tribution. 

The image of the Bullet cluster shown in Fig. [U which 
has unresolved continuum sources subtracted, was con- 
volved to a final beam of FWHM 23'.'8 x 21'.'4 at 80° P. A., 

^ The 1.344 GHz image with embedded unresolved sources sub- 
tracted was kindly provided to us by Dr. Liang. 
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Figure 2. Grey scales represent the 18 GHz mosaie image of 
the buUet cluster with unresolved continuum sources subtracted; 
the image has been smoothed to a beam of FWHM 30". Over- 
laid are contours representing the expected radio halo component 
in the 18-GHz ATCA mosaic image; contours are at 20, 40, 60, 
and 80 /jjy beam ~^. Locations of discrete sources detected by 
ILiang et al I II2000I ) and in this image are marked with cross sym- 
bols. 



which is the beam of the 1.344 GHz image. The halo im- 
age was then subtracted from this image and a 18 GHz 
image of the Buhet cluster, with all continuum sources 
including the halo subtracted, is in Fig. [3] (upper panel). 

4. A DISCUSSION OF THE SZE IN THE BULLET 
CLUSTER 

The sky distribution in the SZE does indeed coin- 
cide with the radio halo distribution in the western and 
southen parts; however, SZE is not detected towards the 
eastern end of the radio halo. Although the depth of the 
SZE detection may suffer errors owing to uncertainties in 
the extrapolation of the halo to 18 GHz, it is clear from 
our Fig. |3] (upper panel) that the detection and extent of 
the SZE is established. Assuming an uncertainty of 0.1 
in the spectral index of the halo, we infer an uncertainty 
of ±20% in the peak SZE dip and 5" in its position. Over 
most of the SZE distribution as well, the estimated error 
in SZE decrement is ±20%. 

We may derive a temperature scale for the SZE image 
using the observing frequency and beam FWHM: the im- 
age rms noise of 8.4 /iJy beam~^ corresponds to 35.3 ^K. 
The deepest SZE feature is at RA: 06'^58"^32!3, DEC: 
-55°57'03", and is clearly offset 35" to the SE of the 
eastern X-ray emission peak; this deep SZE hole is a flux 
density decrement of — 186/iJy/beam and corresponds to 
a brightness temperature decrement of 0.78 niK. The ex- 
pected SZE temperatur e decrement given a C ompton-?/ 
parameter of 3.3x10"'' (IHalverson et al.ll2009D is ATg'z = 
-1.983 y TcMB = -1.79 mK (jHalverson et al] (|200a) es- 
timate —0.77 mK as the peak SZE decrement at 150 GHz 
from their APEX-SZE observations). Our observed SZE 
decrement of —0.78 mK at 18 GHz is significantly lower. 

It may be noted here that our 18 GHz ATCA image 
of the SZE distribution is an interferometer image and 
would miss larger scale extended features. The smallest 
projected baselines in the ATCA observations correspond 
to the antenna diameter of 22 m, and the largest angular 
scale structures that may be reproduced in the ATCA 



image would be less than the antenna beam FWHM of 
2.'6. Because of missing large-scale SZE structure the 
absolute depth of SZE features might be reduced. Sig- 
nificantly differing depths in the SZE decrement in the 
APEX and ATCA images imply that there is substantial 
extended SZE structure in the Bullet cluster on angular 
scales of 2.'6 and larger, which is resolved out and missing 
in the ATCA inteferometer image. 

In Fig. [3] (bottom panel) we show contours of the SZE 
overlaid on Chandra X-ray intensity distribution shown 
using grey scales. 
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Figure 3. Upper Panel: Grey scales and contours represent the 
18 GHz mosaic image of the SZE in the Bullet cluster; the image 
has been smoothed to a beam of FWHM 30". Contours are at 
(—16, —12, —8, —6, —4, —3, —2, 2, 3) times the image rms noise 
of 8.4 /xjy beam~^. The peaks in ROSAT HRI X-ray emission are 
marked with cross symbols. 

Bottom Panel: Contours represent the 18 GHz mosaic image of the 
SZE in the Bullet cluster; the image has been smoothed to a beam 
of FWHM 30". Contours are at (-16, -12, -8, -6, -4, -3, -2, 
2, 3) times the image rms noise of 8.4 /^Jy beam~^. Grey scales 
represent Chandra X-ray emission. 



1. We observe a deep SZE 'hole' that is significantly 
displaced to the SE of the X-ray peak in the hot- 
ter eastern cluster. Since the SZE probes the inte- 
grated pressure along the line of sight, this suggests 
that there are significant pressure variations in the 
intracluster gas, presumably because of the ongoing 
cluster collision. X-ray spectral luminosity is pro- 
portional to / nlT~^dl, whereas the SZE is propor- 
tional to / neTdl: the two effects represent different 



weightings of temperature T and electron density 
no. X-ray emission images are more representa- 
tive of rif. distribution, whereas the SZE images are 
representative of T distribution. The offset of the 
SZE 'hole' from the emission peak clearly suggests 
that the peak in intracluster gas temperature is 
spatially offset from the peak in gas density in the 
Bullet clustei0. Our observations are the first time 
such an effect has been observed; the high angular 
resolution ATCA observations of the Bullet cluster 
presented herein may be the first time that such 
a detailed comparison between the SZE substruc- 
ture and X-ray emission has been possible in any 
cluster. The offset is a demonstration of the com- 
plexity in intracluster gas distribution properties in 
merging clusters. Apart from the deep SZE 'hole', 
we observe significant SZE substructure across the 
cluster complex. All these SZE components are 
displaced from the peaks in X-ray emission corre- 
sponding to the two clusters comprising the Bullet 
cluster complex. 

2. All substructure we detect in the SZE image lies 
within the boundary of the radio halo. Moreover, 
careful examination of Fig. [2] indicates that the SZE 
in the western parts of the cluster complex appear 
to follow the contours of the radio halo. This is sur- 
prising, since the radio halo arises from synchrotron 
plasma and the SZE from thermal intracluster gas 
and there is no spatial correspondence expected be- 
tween these distinct components. A plausible rea- 
son for the cospatial locations of the SZE and ra- 
dio halo — in the western parts of the cluster — is 
a common origin for the two populations of ra- 
diative electrons. Radio halos are believed to be 
synchrotron emission from electrons accelerated (or 
re-accelerated) to relativistic energies in the turbu- 
lence generated by cluster mergers. It may be that 
the radio halo was created in the turbulent wake of 
the passage of the 'Bullet' through the cluster gas, 
and the SZE also arises cospatially from relatively 
hotter electrons in the same turbulent wake. Ab- 
sence of SZE towards the eastern parts of the radio 
halo would, in this picture, indicate longer lifetimes 
for the rela tivistic electrons compa red to the SZE 
electrons. (jMarkevitch et ahl 120021 ) have reported 
finding spatial correlation between halo brightness 
and local gas temperature in merging clusters; this 
is supplementary to the SZE and the halo being 
cospatial in our observations. 

3. The SZE decrement is relatively shallow towards 



the peaks in the X-ray emission. In Fig. |3] it may 
be noted that the SZE contours are indented in sky 
regions where the X-ray emission is most intense. 
The observed avoidance of the highest emission 
measure regions by SZE suggests that the high- 
est density parts of the cluster gas may have some- 
what cooler temperatures, owing to shorter cooling 
timescales, and consequently contribute relatively 
less to the SZE. 
Wc note here that Massardi et al. (ApJL submitted) 
found a similar displacement of SZE dip from X-ray 
brightness peaks in the merging cluster CL0I52— 1357. 

5. SUMMARY 

We have presented 18 GHz observations of the Bullet 
cluster from the ATCA, and report robust detections of 
substructure in the Sunyaev Zeldovich Effect. We sub- 
tract unresolved continuum sources in the cluster field 
using 18 GHz ATCA images made with enhanced weight- 
ings to the longer visibility spacings. Wc make reason- 
able assumptions about the radio halo in the cluster to 
extrapolate its observed properties from 0.8-8 GHz to 
18 GHz and model and subtract the radio halo contribu- 
tion. We image substantial substructure in SZE towards 
the Bullet cluster, finding a deep SZE hole that is signif- 
icantly displaced from the peak in X-ray emission. The 
main conclusion from this work is that the intracluster 
gas in the merging Bullet cluster appears to have signifi- 
cant pressure distribution structure and peaks that differ 
from that in gas emission measure, galaxy and dark mat- 
ter distributions. The work emphasizes the complexity 
in merger dynamics, and suggests that modeling the SZE 
contributions from cosmological clusters — for SZE clus- 
ter counts as well as SZE contributions to small-scale 
CMB anisotropy — is nontrivial. 

The Australia Telescope Compact Array is part of the 
Australia Telescope which is funded by the Common- 
wealth of Australia for operation as a National Facility 
managed by CSIRO. 
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